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ABSTRACT 



Global temperatures are obtained using a linear least 
squares regression method with satellite radiation measure 
ments, particularly those obtained from the SIRS-B aboard 
the NIMBUS IV satellite. Regression equations relating 
temperature to spectral radiance observations are employed 
The regression equations were then applied to independent 
observations to determine the feasibility of temperature 
determinations over sparse data regions. 
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I. 



INTRODUCTION 



The launch of the NIMBUS IV on 8 April 1970 was the 
second step in a process of developing a means of determin- 
ing the three-dimensional state of the atmosphere by indir- 
ect means. The orbit of the NIMBUS IV was a near-circular, 
polar-orbit, with a height of approximately 600 nm . The 
new instrument, known as SIRS-B (which is similar to the 
instrument on NIMBUS III, known as SIRS-A)) has the added 
capabilities of "side looking" scan for greater area cover- 
age, and the determination of the tropospheric water vapor 
content. 

The calibration procedure used to transform the radia- 
tion sensed in the eight channels of the SIRS-B to spectral 
radiances is discussed in the NIMBUS IV USERS GUIDE (1970). 

The temperature profile of the atmosphere within the field 

2 

of view, comprising an area of (225km) at the subsatellite 
point, is deduced using the channel spectral-radiance en- 
semble. The particular wavenumbers sensed are small inter- 
vals in the atmospheric window and the 15-micron band of 
C02* The channel numbers and wavenumbers are listed: 



CH. .NO. 
1 


1 2 


3 


4 


5 


6 


7 


8 


WAVE. NO. 


899.0 750.0 


734 . 0 


709 . 0 


701.0 


692 . 0 


679 . 8 


668.7 
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Through the equation of radiative transfer, the meas-- 
ured radiances are related to the temperature-structure 
T(p) of the atmosphere, the temperature at the earth’s sur- 
face and that at the cloud tops which may be within the 
f ield-of -view . Assuming that a single cloud layer exists, 
the following set of integral equations can be constructed 

for the measured radiance N(v.) in channel i: 

T 

N(V^)=A j B[v^,T(p^) ]t(V^,P^) ‘^BEv^.TCp) ]dx(V^,p) ) 

( 1 ) 

S 

+ (l-A) j B[v^T(p^)] T(v^,p^) -j/ B[v^,T(p)]dT(v^,p) ) 

i=l, 2, , . . , 8 is the channel number, 

N(v^) = spectral radiance at wavenumber 
B[v^,T(p)] = Planck radiance at wavenumber and tempera- 
ture, (Eq . (3) below), 

t(v.,p)= fractional transmittance of the atmosphere 

in the spectral interval centered at wave- 

number V. and assumed known as a function of 
1 

the pressure p at the radiating level, 

A = the product of the fraction of cloud cover 

within the f ield-of -view and the cloud emis- 

sivity (which is usually assumed to be unity). 

Subscripts c and s refer to cloud-top and earth surface. 

The special case of equation (1) where A=0 gives the 

dear-column or corrected radiance N (v . ) , This term is 

c X 

contained within the second brace of equation (1), It is 
this quantity which must be determined in all channels in 
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order to convert the corrected radiance set into a vertical 
temperature profile T(p). The clear^column spectral radi- 
ance at the top of the atmosphere at wavenumber is there- 
fore : 



Nc(Vi) = BIv^,T(p^)]t(v^,p^) 



where 



B[Vj,,T(p)] dT(v^.pg) 



B[v^,T(p)] = {exp [c 2 V^/T(p)] 




( 2 ) 



(3) 



is the Planck function at wavenumber . 

Here T=T(p) is the temperature (assumed to be in thermo- 
dynamic equilibrium at level p) , and c^ and are known 

cons tan ts 

-5 2 -1 -1 

c^ = 1.9061 X 10 erg cm sec ster 

C 2 = 1 . 43868 cm 

x(v^,p) is the transmittance at wavenumber of 

the atmosphere above level p. 

Smith, Woolf and Jacobs (1970) have used equations (2) 

and (3) to de f ine a weighted-mean black-body temperature 

. in terms of the corrected radiance N (v . ) which would 
B , 1 C l 

be observed from an atmosphere of structure T=T(p) 






(4) 



Here T_ . is independent of pressure. Solving the relation- 

D 1 



ships (3) and (4) for T . gives 

D 1 



T (v ) = c V /In j Ic V ^/N (v.)] + 1 j 

D X zi iJLlci ^ 



(5) 



Smith et al (1970) have provided multiple regression equa- 
tions which relate NMC effective grid point temperature 
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T(p.) at SIRS-A scan spots to the sixteen predictors T (v . ) 

J B r 

2 

and Tg(v^) • They have done this by relating a two*-week 
NMC data file of temperatures at standard pressure levels 
p^ to the corrected radiances, and arrived at an equation 
of the form 

8 _ 

T(p ) = T(p ) + Z a(v. ,p.) [T_ (V.)- T^(V.)] 

J J ^ ^ B r 15 1 

® — 2 
+ 2 a’(v.,p.) [T^(V.)- T_(V.)]^ (6) 

i=l ^ J B 1 B 1 

A similar study will be made for T(500) and T(300) using 
a data-base of corrected radiances from SIRS-B and NESC de- 
rived temperature profiles T(p^) together with simultaneous 
FNWC ef f ec tive grid-point soundings (available from the 
Fleet Numerical Weather Central ADP tapes) at locations of 
the scan-spots in both space and time. For the SIRS-B data 
file, a revised NESC formula similar to equation (6) was 
available for the conversion of corrected radiance-sets 
into a statistically consistent profile T=T(p). 
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II. DATA PROCESSING 



A. THE ORIGINAL DATA 

The original data utilized in this paper were provided 
by the National Environmental Satellite Center (NESC) , and 
consisted of computer listings of a selected set of SIRS-B 
scan spots recorded during December 1970. In addition to 
the geographical coordinates and time of each observation, 
there were also listed the raw (un c or r e c t ed ) radiances in 
each channel. The raw radiances are assumed to be derived 
from equation (1) for the one cloud-layer model and the 
corrected radiances resulted from the systematic elimina- 
tion of cloud contamination effects from each channel [fol- 
lowing a technique of Smith et al (1970)]. 

The atmospheric window channel (channel 1) was inopera- 
tive during this period; therefore, the values of the cor- 
rected window channel were chosen from locations where the 
surface temperature, based on the NMC data file, was known 
to be at a pressure-level close to p = 1000 mb, where T^ is 
known. Furthermore, since the radiance-correction computa- 
tion involves iterative use of equations (2) and (6) (using 
an equivalent ”clear-column*’ temperature profile, T = T(p), 
at each iteration) , all original radiance-sets chosen could 
thus be "tied down" by the selection of scan-spots having 
surface pressure values p^ = 1000 mb, with T^ known. 



/ 
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It should be noted that the data samples provided were 
enhanced by the side-viewing capability of the SIRS-B. The 
side-viewing capability enabled the scanning of more atmos- 
phere-earth columns both on and off the subtrack for which 
the approximation p^ = 1000 mb was valid. Slant-columns 
thus viewed are useful for reduction to dear-column tem- 
perature profiles T = T(p) after the slant path is normalized 
to a single atmospheric depth over the spot viewed. 

In general, virtually all scan spots have some effective 

2 

cloud cover (A > 0) in view of the large area (225 km) be- 
ing viewed. Of the seven CO 2 channels being used, the two 
channels subject to the maximum attenuation by the inter- 
position of a "cool!* cloud top within the desired T = T(p) 
profile are channels 2 and 3. 

Smith, Woolf, and Jacobs (1970) have indicated the pro- 
cedure involved for correcting for cloud-contamination of 
the eight channels of radiances. Based upon the actual NMC 
climatological update of each of the standard pressure- 
levels, (10, 30, 50, 100, 150, 200, 250, 300, 400, 500, 

700, 850, and 1000 mb) they have derived multiple regression 

equations relating T(p^) to the clear-column-radiance 
N^(v^), or to the blackbody temperature-equivalent of the 
dear-column radiance. For example. Smith et al (1970) 
obtained regression equations essentially of the form of 
equation (6) and have listed the least-squares best-fit 
coefficients obtained for a period in September 1969 appli- 
cable in a zonal band 35-55°N. As a second approximation, 



14 



Smith et al have also generated a six-predictor analog of 
equation (6) with channels 2 and 3 deleted for the same 
time frame as utilized for the eight-predictor equation. 

In both cases, they have deduced equation (6) and its six- 
predictor analog by employing a large sample-size of simul- 
taneous values of corrected radiance-sets and of values of 
T(p^) (at standard pressure levels), interpolated from the 
NMC file to the location of the SIRS scan spot. The six- 
predictor and eight-predictor regression equations are con- 
tinually updated to a revised two-week data base after 
three additional data-days are accumulated. Such equations 
were also available during the SIRS-B period of December 
1970. 

The procedure for correction of the ^th channel radi- 
ance makes iterative use of the updated SIRS-B operational 
six-predictor regression equation in deducing a first guess 
of T=T(p). At the first guess, the uppermost five radiances 
(N^,...,Ng) which are least affected by cloud contamination 
are used as the first guess correct-radiance set, together 
with the prescribed from the window-channel blackbody 

*'fix'* at the scan spot. With these six-predictors known, a 
first guess of the temperature profile, T(p^), at all 
standard levels is computed. After the first-guess of the 
temperature profile is made, the profile is then applied to 
channels 2 and 3 to determine the **best fit” for cloud 
height and effective cloud cover A using equation (1), re- 
written for simplicity as 
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(7) 



N(V.) = (1--A) N^[v^,T(p) ,p^] + AY [V.,T(p),p] 

Here Y represents the contents of the first brace of equa- 
tion (1). Also gCv^) is defined as 

c(Vi) = (V.) - N(V.) (8) 

for the correction to be applied to the radiances for chan- 
nels 4 through 8, Writing equation (7) for channels 2 and 
3 in terms of the corrected values of based on T(p) (the 
dear-column temperature) and the corresponding computation 
of Y from the top of the clouds gives two simultaneous 
equations. From these two simultaneous equations, one 
solves for the two unknowns, A and p^, using indirect 
methods for solving for p^ (i.e. making separate tests for 
p^ = 200, 250, 850 mb) and, in turn, determining the 

value of A which ”best-fits" equation (7) expressed in 
terms of N(V 2 ) and NCv^)* These values of A and p^ are 
substituted into equation (7) for a second corrected esti- 
mate of N, , . . . . ,N^and a new temperature profile is calcu- 

H O 

lated. The iterative procedure is continued in this manner 
and usually converges within four iterations, i.e., the 
correction c(v^) of (8) usually converges 
cCVf) liinit) 

for a suitable tolerance within four iterations. At this 
point, the finalized temperature profile T(p) may be used 
to compute the dear-column radiances for channels 2 and 3 
as well as the five uppermost CO^ channels using Eq. (I)* 
Cloud influence is thus compensated when the final choice 
of temperature profile has been made. 
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In a sizeable fraction of the cases selected, the con- 
vergence c(v^) required more than four iterations which 
means that a relaxed requirement for convergence had to be 
employed. This generally occurred with cloud tops at or 
above 500 mb. Such cases of below-normal attenuation re- 
covery were indicated on the NESC da ta-lis tings as decreased 
reliability. Nevertheless, no data cases were discarded 
for such reasons. 

B. ADJOINING OF FNWC DATA 

Analyzed temperature-field data for 500 and 300 mb cover- 
ing the period from 00 GMT 24 December 1970 through 00 GMT 
29 December 1970 was provided by Fleet Numerical Weather 
Central (FNWC). The FNWC temperature data had not been 
correlated with radiances as was the NESC data. The FNWC 
temperature fields at 500 and 300 mb [T (5) and T (3) ] were 

r r 

spatially and temporally interpolated to obtain consistency 
with the SIRS soundings. Spatial interpolation was accom- 
plished using Bessel’s central-difference interpolation 
[Haltiner, 1971] to the latitude and longitude of the SIRS-B 
soundings. Temporal interpolation to the time of the scan- 
spot was accomplished utilizing linear interpolation between 
temperature fields preceding and succeeding the time of the 
SIRS sounding. These methods were employed to obtain 

statistical dependence of T„(5) and T^(3) upon the corrected 

r r 

radiance set N-,....,Nq. This necessitated the printout of 
1 o 
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FNWC interpolated temperatures in a format consistent with 
the requirements of the program BIMED 02R, to be used in 
(III). 

C. DATA DESIGNATION 

The original NESC data comprised some 207 to 329 scan- 
spots .per day and were identified by the specific time of 
the SIRS soundings and their latitude and longitude. The 
original FNWC data tapes consisted of temperatures at the 
levels 300 and 500 mb, at each grid point of the 63 x 63 
Northern Hemisphere FNWC grid. The FNWC data at scan spots 
were obtained by interpolation as previously described. 

A three-day data base of both NESC and FNWC data was 
thencompiled using data from 0600 GMT to 1800 GMT for 24, 
25, and 26 December 1970. A four-day data base was com- 
piled by the addition of data from 0600 to 1800 GMT 27 
December 1970 to the three-day data base. Data sets for 
independent regression tests were compiled using data cen- 
tered on 00 and 12 GMT 27 December and 00 and 12 GMT, 28 
December (i,e,, data-base times plus 12 and 24 hours). 

All data sets were divided into three bands. These 
bands were bounded by latitudes: 20-40°N, band 1; 35-55°N, 
band 2; and 50-70°N, band 3, The bands were overlapped to 
reduce the possibilities of discontinuities across the 
boundaries when regression calculations were made. 
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The designation of dependent and independent variables 
was made and the dependent variable was designated T(i,j) 
and the independent variables were designated as shown; 



CH. NO. 12345678 

WAVE NO. 899.0 750.0 734.0 709.0 701.0 692.0 679.8 668.7 

VAR. NO. N, N. N- N, N^ N, N., N. 

12345678 

For each scan spot selected in the test period, the 

time and location were encoded on IBM cards, as were the 

temperatures T„(5), T (3), T (5), T (3) . Here, the sub- 
N M r r 

script N indicates an NESC-genera t ed temperature, using the 
convergent procedure described previously, and the subscript 
F denotes interpolated FNWC temperatures. Also encoded 
onto IBM cards were the eight channels of corrected radi- 
ance s ( provided by NESC using the method described in II, A.) 
for each of the scan spots. All data were encoded in a 
manner consistent with the requirements of the stepwise 
regression program, BIMED 02R (Dixon 1966) which will be 
discussed later. 
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Ill, THE REGRESSION HETHOD 



A, THE STEPWISE REGRESSION PROCEDURE 



Many problems in research require the extensive analysis 
of large amounts of data. The data handling process should 
be made as automatic and rapid as possible. The appropriate 
tools for analyzing large data samples with numerous inde- 
pendent variables often require multivariate statistical re- 
gression techniques to attack the problem. The Biomedical 
Computer Programs (Dixon 1966) were developed to provide 
the techniques necessary for the statistical and mathemati- 
cal analyses required. One of these programs is the Step- 
wise Regression Analysis Program BIMED 02R. 

BIMED 02R computes, in a stepwise manner, a sequence of 
multiple linear regression equations. One variable is added 
to the regression equation at each step. The variable added 
is that which accomplishes the greatest reduction in the 
previously unexplained sum of squares. This variable is 
also the one which has the highest partial correlation with 
the dependent variable at the particular step in the analy- 
sis of the variance. Equivalently it is the variable which 
would have the highest F-value if it were added. The F- 
value, or F-statistic upon entry, F^^ is expressed at step k 
as (Dixon 1966) 



F^(l,n-k-l) 



%(C.E.V.k) - %(C.E.k-l ) 
%(U.E.V. ,k) 
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where 



%(C.E.V.k) is the percent cumulative explained variance, 
step k. 

%(C.E.V.k-l) is the percent cumulative explained 
variance, step k-1. 

%(U.E.V.k) is the percent unexplained variance remain- 
ing at step k • 

This study utilized a statistical model expressible in 
the form 

T = C H-C,N,+C.N. + C.N. + C,N,+C^N^ + C,N, + C^N^ + C^Nq (9) 

01122334455667788 

where C through are to be determined by the stepwise 
o o 

least squares technique. Equation (9) may be expressed in 

the more compact form utilizing the matrix-product equation 

1 

Ni 

T = (c c, , c ) ( : ) (10) 

8 

In (10), the notation T^(3,5) will be used when the NESC 

data at 300 or 500 mb (respectively) are correlated against 

the predictors (N- >Nq); whereas T_(3,5) is the nota- 

i o r 

tion for the corresponding predictands with FNWC tempera- 
tures. Equation (10) may also be used in application with 
independent data occurring up to 12 or 24 hours later than 
the dependent sample which yielded the particular coeffi- 
cient matrix (C , ...., C^). 

o 1 o 

B. RELATED STATISTICAL PARAMETERS 

Available as outputs from the BIMED 02R program are 
"several related statistical parameters. Among these sta- 
tistical parameters are: 
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a. multiple R 

b, standard error of estimate, S. E. 

c. mean value, T 

d, standard deviation a 
e • F-value 




The significance of R in relationship to S, E, is as 
follows (Crow et al 1955) : 

(S.E.)^ = 0^ [ (n-l)/(n-k-l) ] (1-R^) (11) 

where 

n _ ^ 

= Z (T. - T) /n 
i = l ^ 

is the variance of the temperature sample, and 
n = sample size 
i = sample-element index 
k = number of predictors 

Furthermore, R is the multiple correlation coefficient after 

k predictors are added, with emphasis on k=8, in order to 

derive maximum information from the SIRS-sounding . It fol- 

2 

lows that the percentage unexplained variance (1-R ) for 
sample sizes in the range 171 to 228 as contained in the 
dependent data sample are very closely approximated by 
(S.E./a)^ = 1-R^ 

2 

The percentage explained variance is then given by R ; that 
is 

R^ = 1 - (S. E./a)^ (12) 

s^ince (n-k-1) / (n-1) is close to unity for the dependent 
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sample sizes considered here* The factor (n-k-1) / (n-1) 
becomes (n-2)/(n-l) = 1 for the independent test cases 

considered since in this case, there is only one way of 
forming the predictor to be listed* 

C. THE REGRESSION PROCEDURE APPLIED TO THE DEPENDENT 
SAMPLES 

The test on the dependent data was made using the least- 

squares, stepwise regression technique of BMD 02R to develop 

the coefficients C ,C^,..*.,Co* The coefficients for the 

o 1 o 

three-day and four-day dependent samples for the T(5) and 

T(3) regression equations were computed by use of this 

technique. This test was performed on both the NESC-regres- 

sion generated data and the interpolated FNWC data* The 

coefficients for the NESC three-day data sample are listed 

in Table 1, part (a) for 500 mb; and part (b) for 300 mb* 

Table 2 contains similar results for the four-day sample 

base of NESC data. The coefficients for the FNWC data are 

listed in Table 3 (three-day case) and Table A(four-day). 

The relevant statistics for the tests on the several 

dependent data bases are listed in Table 5 and 6. The 

statistics included in these tables are the mean, standard 
2 

deviation, R , and standard error. A discussion of the 
statistics of prime interest is provided in IV. 
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TABLE 1. Regression coefficients for NESC dependent data. 

(a) 500 mb, three-day sample 



CONST 


20-40N 


35-55N 


50-70N 




-43.1058 


-52.2372 


-90.9234 




-0.0327 


0.4187 


0.1508 




-0.0932 


-0.4795 


-0 . 2562 




0.8922 


-0.2670 


0.0 


^4 


0.6949 


3.0232 


2.9128 




-1 . 8871 


-2.6461 


-1 .7719 


"6 


-0.9188 


-0.0726 


-0.1134 




0.4129 


0.1608 


0 . 2351 




0.6664 


0 . 0764 


-0.1635 




(b) 300 mb, 


three-day sample 


o 

o 


-70.5523 


-113.6468 


-108.4219 




0.4109 


0.5032 


0.3409 


S 


-0.2953 


-0.3192 


-0.2337 




-1.4766 


-1.4030 


-0.6563 


^4 


3.0286 


3.1138 


1.3928 




-0.7404 


0.0857 

* 


1.0938 


"6 


-0.5295 


o 

o 


0.1261 




-0.5835 


-1.0623 


-1.1938 




0.6865 


0.4809 


0.3922 




indicates that this 


variable (*) 


was not accepted 




because it failed to 


explain appreciably 
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TABLE 2. Regression coefficients for NESC dependent data. 





(a) 500 


jnb , 


four-day s 


ample 


CONST 


20-40N 




35-55N 


50-70N 


S 


-48.6906 




-51.8292 


-92.4865 




0.0 




0.4842 


0.1143 


S 


-0.0907 




-0.6454 


-0.1644 


s 


0.8032 




-0.1472 


-0.0980 




0.6590 




3.0011 


3.0117 


s 


-1.5652 




-2.6431 

* 


-1 .7844 




-0.9017 




0 

• 

0 


-0.1327 


"7 


0 .2990 




0.1605 


0.1894 




0.6797 




0.0354 


-0.1311 




(b) 300 


mb , 


four-day s 


ample 




-101. 7840 




-108.5859 


-108.1667 


"1 


0.1288 




0.5599 


0.3757 




-0.0794 




-0 . 5138 


-0.3326 




-0.2788 




-1 . 0924 


-0.5185 




0.5219 




2.7352 


1.2617 


s 


1.2700 




0.1432 

* 


1.1646 


s 


-1.0121 




0 . 0 


0 .1026 


s 


-0.5885 




-1 . 0406 


-1.2316 


s 


1.0295 




0.5160 


0.4386 




* 

indicates that 


this 


variable 


(*) was not accepted 



because it failed to explain appreciably 
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TABLE 3. Regression coefficients for FNWC dependent data 



(a) 500 mb, three-day sample 



CONST 


20-A0N 


35-55N 


50-70N 


S 


-15.6121 


-55.3823 


-85.1272 




0.3605 


0.4119 


9.0638 


s 


-0. 8342 


-0.8714 


-0.7219 


s 


0.7760 


0.9170 


1.7073 


"a 


1.2602 


1.6287 


0.6259 


s 


-1.8509 


-1.9698 


-0.6344 


"6 


-0.9308 


0.1453 


-0.1520 




0.2310 


-0.1858 


0.0* 


S 


0.4440 


0.1834 


-0.0545 




(b) 300 


mb, three-day sample 




S 


-56.7715 


-107 . 1405 


-99.6836 




0.4459 


0.2963 


0.3651 


S 


-0.6807 


0.0459 


-0.6303 


s 


-0 . 3890 


-1.1144 


0.3020 




1.9185 


2.1166 


0.5930 




-0.5798 


0 . 3357 


0.9944 


^6 


-0.9211 


-0.1590 


0.0444 




-0.2675 


-0.5875 


-0.6276 


"8 


0.5916 


0.2311 


0.0392 


* 

indicates 


that this variable (*) was not 


accepted 


because it 


failed to 


explain appreciably 
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TABLE 4. Regression coefficients for FNWC dependent data. 

(a) 500 mb, four-day sample 



CONST 


20-40N 


35-55N 


50-70N 




-15.3487 


-51.3814 


-88.3277 


"l 


0.4477 


0.5272 


0.0260 




-1 .0 928 


-1.2317 


-0.5774 


s 


1.0306 


1 . 3415 


1.4995 


"4 


1.1163 


1.3383 


0.8176 




-1 .7344 


-1.9340 


-0.6736 


^6 


-0.8981 


0.2096 


-0.1147 




0 . 2645 


-0.1328 


-0 . 0979 


^8 


0.3767 


0.1141 


^ -k 

0 . 0 




(b) 300 mb 


, four-day sample 






-53.5409 


-103 .0841 


-99.9952 




0.4850 


0.3744 


0 . 3821 


S 


-0 .7624 


-0.1827 


-0.6721 


S 


-0.4110 


-0.8859 


0.3722 


^4 


2 .0702 


2.0795 


0.5005 


S 


-0.6480 


0.1929 


1.0677 


^6 


-0 . 8955 


-0.1818 


0.0171 


S 


-0 . 1607 


-0.5045 


-0.6575 




0 . 4292 


0.2066 


0.0757 


* 


indicates that this 


variable (*) was 


not accepted 



because it failed to explain appreciably 
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IV. DEPENDENT DATA RESULTS 



The detailed results of the tests on dependent data are 
available in Tables 5 and 6. However, the essential por- 
tions needed for comparisons by the different tests using 
dependent data have been listed in Table 7 where the values 
of fractional explained variance as originally listed in 

Tables 5 and 6 are summarized. Also listed in Table 7 are 

2 

the values of R averaged over the data-base specification 

2 

period (or mode) for each band as well as the values of R 
averaged over the bands for a given data-base mode. Com- 
parisons between regressions were based on the fractional 

2 

explained variance R rather than on the values of standard 
2 

errors, as R has been normalized with respect to the sam- 
ple variances, which vary between bands. 

The third value listed in the band columns of both parts 
(a) and (b) of Table 7 is the average over the data-base 
modes in each band. The values listed in the mode-average 

row and band-average column are the overall averages by the 

2 

statistical regression. Comparisons of R are made of the 

overall averages and of band averages, as well as data-base 

mode averages at both the 500-and 300-mb levels. Differ- 

2 

ences in values of certain compared R are termed ”% shrink- 
age** in the explained variance arising from the specifica- 
tion tests under comparison. 
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A. COMPARISON OF NESC SPECIFICATION STATISTICS 
Examination was first made of the left hand side of 

Table 7 for comparison of the NESC regression results at 

300 mb relative to those at 500 mb. The outstanding 

2 

result of this comparison reveals R^(3) to be smaller than 

2 

R^(5) by 0.1660, or 16.6% when averaged over bands 1, 2, 

and 3. Upon comparison of individual band results (averaged 

over the three and four day specification modes) one finds 

the same order of shrinkages distributed over the bands, 

namely, an average of 16.6%, but with a larger than normal 

2 

shrinkage in band 2 where the value of R^ fell off to 
0.6857 . 

Upon use of the four-day estimation equation (10) in 

the respective cases (left side of Table 7) there was a 

small shrinkage in explained variance (relative to the 

three-day results) at both the 300- and 500-mb levels. 

This result was not quite anticipated since Smith et a] 

found it advisable to use data-banks of two weeks for 

specification. Nevertheless, the small differences between 
2 

the R^ for the two prediction-base methods, whether at 300 
or 500 mb, were not large enough to be considered signifi- 
cant in the tests conducted here. 

B. COMPARISON OF FNWC SPECIFICATION STATISTICS 

Here our comparisons will mainly concern the right hand 
side of Table 7. Again, predictor-mode differences seem to 
be marginal so that for the present they will be disregarded. 
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2 2 

Comparison of Rp(3) with Rp(5) shows that band average 

2 2 

shrinkage in the amount R„(5) - R„C3) = 0.1480 occurs. A 

r Jb 

breakdown of this statistic by bands shows shrinkage in the 
amounts 9.16%, 22.39%, and 12.87% for bands 1, 2, and 3, 
respectively . 

Comparisons made between NESC and FNWC regression re- 
sults (left side and right side of Table 7) show approxi- 
mately 2% difference between R^(5) = 0.8737 and R^(5) = 

r In 

0.8957 and less than 0.5% difference at the 300-mb level. 

Thus, the FNWC specifications at both levels explained very 

nearly as much variance as did the NESC results. There 

were some variations between bands, but generally, the 

same geographic pattern persisted throughout: namely, 

2 

moderate to high R in bands 1 and 3, with the lowest value 
2 

in R occurring in band 2. This was true of all FNWC speci- 
fications (300 and 500 mb) as well as of the NESC 300 mb 

specification. The NESC 500 mb specification gave rise to 

2 

an increase in R with an increase in latitude. 

In summary, both NESC and FNWC regression equations 
tested well in terms of explained variance at 500 mb. Both 

showed diminished specification at 300 mb to values of 

2 . 

R = 0.73. Finally, in no case was there a significant 
change when testing the four-day data base results against 
those for the three-day data base. 
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TABLE 5. Dependent-data temperature specification 
statistics resulting from pooled samples of NESC data 



Band 


(a) 

Sample 

Size 


500 mb, three-day 
Std, 

Mean Dev* 


sample base 


Std. 

Error 


20-A0N 


171 


-15.5029 


6.8730 


0.8533 


2.6970 


35-55N 


191 


-28.3298 


6.9763 


0.9113 


2.1218 


50-70N 


177 


-3A . 1A69 


7 . 0397 


0 .9321 


1. 872A 




(b) 


500 mb, four-day 


sample base 




20-A0N 


206 


-15 . A369 


6 .7276 


0.8379 


2 .2275 


35-55N 


228 


-28.1623 


6.7376 


0 .90A0 


2 . 1202 


50-70N 


216 


-3A . 2176 


7 .0982 


0.9358 


1.8336 




(c) 


300 mb, three-day 


sample base 




20-A0N 


171 


-A2 . 7895 


A .6668 


0.7719 


2 . 283A 


35-55N 


191 


-50 . A869 


A . 7A27 


0.68A2 


2.7152 


50-70N 


177 


-5A .6667 


3.9610 


0.7833 


1.887A 




(d) 


300 mb, four-day 


sample base 




20-A0N 


206 


-A2 .7233 


A .6398 


0.65A8 


2.7809 


35-55N 


228 


-50.3509 


A . 6875 


0.6871 


2 .6637 


50-70N 


216 


-5A.7778 


3 .9610 


0.7970 


1 . 8186 
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TABLE 6. Dependent-data temperature specification 
statistics resulting from pooled samples of FNWC data 





(a) 


500 mb, three-day 


sample base 






Sample 




Std . 




Std . 


Band 


Size 


Mean 


Dev . 




Error 


20-40N 


171 


-15.9520 


6 .6156 


0.8760 


2.3861 


35-55N 


191 


-28.2174 


7.3512 


0 . 8704 


2 . 7045 


50-70N 


177 


-34 . 2333 


7.3384 


0 . 8816 


2.5767 




(b) 


500 mb, four-day : 


sample base 




20-40N 


206 


-15.9179 


6.4674 


0.8651 


2.4232 


35-55N 


228 


-28.0130 


7 . 1574 


0.8598 


2 .7283 


50-70N 


216 


-34.2405 


7 .3470 


0.8894 


2.4842 




(c) 


300 mb, three-day 


sample base 




20-40N 


171 


-42 . 2974 


5.2983 


0.7833 


2.5267 


35-55N 


191 


-50 . 8692 


4.6688 


0.6413 


2.8571 


50- 7 ON 


177 


-54.9912 


4.3548 


0.7494 


2.2313 




(d) 


300 mb, four-day 


sample base 




20-40N 


206 


-42 . 3044 


5.1558 


0.7747 


2.4962 


35-55N 


228 


-50 .7638 


4.6363 


0.6411 


2 . 8276 


50-70N 


216 


-55.1295 


4.3142 


0.7641 


2.1354 
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7795 0.7966 0.8620 0.8174 0.8248 0.753 



V, ■ INDEPENDENT DATA TEST RESULTS 



Verification of the dependent regression equations was 
accomplished using the regression coefficients (described 
in III) and the 12-and 24-hour independent data samples. 

The equations of form (10) developed from both the three-day 
and four-day samples were tested. The independent tests 
generate equations of the general form 

T(5,3) = + B^T(5,3) (13) 

where and B^ are regression constants and T represents 

the temperature-estimators from the dependent regression 
equation (10) wi th .coef f icien t s appropriate to the level, 
the band and the mode - s t r a t i f i ca t i ons (Tables 1,2, 3, 4). 

The regression analysis for these tests is similar to the 
techniques applied in IV in that T(5,3) is compared by re- 
gression methods to the observed temperature at each scan- 
spot. The results of the independent tests are given in 
Tables 8, 9, 10, and 11. Table 12 is constructed in a man- 
ner similar to that used to develop Table 7. The major 
difference between Tables 12 and 7 is that the 12- and 24- 
hour independent samples are statistically pooled in appli- 
cation of the independent regression test in order to obtain 
a larger sample size for each band, and also for ease of 
comparisons. The results of the pooled independent data 
verifications are presented in Table 12. 
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Internal comparisons of the values of fractional ex- 
plained variance for the independent data yielded results 
quite similar to those obtained from the dependent data. 

The comparisons are stratified in two manners; by bands and 
by modes, as well as by levels. 

A. COMPARISONS OF NESC INDEPENDENT VERIFICATIONS 

From the left side of Table 12, it is easily observed 

2 

that the largest overall average R = 0.8392 occurred in 
connection with the 500-mb verification of T^(5). The cor- 
responding mean NESC result at 300 mb verified with a per- 

2 

centage explained variance R^(3) = 0.7533, i.e., with a 

2 

shrinkage of 8.59% relative to R^(5). 

As noted in IV (A and B), both NESC dependent specifica- 

2 

tions (300 and 500 mb) in band 3 verified highest with R 

ranging between 0.82 and 0.91. The previously noted trend 
2 

in R continues at the 500-mb level of independent data with 
2 

R increasing poleward from 0.78 in band 1 to 0.91 in band 

2 

3. At 300 mb, this band distribution of R^(3) is not so 
apparent since values close to 0.72 are associated with 
bands 1 and 2 with a mean value of R^(3) - 0.&2 existing in 
band 3 . 

B. COMPARISONS OF FNWC INDEPENDENT VERIFICATIONS 
Comparisons of the FNWC independent tests (Table 12, 

right side) bring out primarily the same results which have 

been developed from dependent base samples. For example, 

. 2 

the FNWC overall R (5) - 0.78 may be compared with 
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2 -V 2 

Rp(3) ^ 0.75, In general, the independent Rp(3,indep) is 

2 

close to the values found for R (3,dep), In keeping with 

X 

earlier results, band 2 has a relative minimum value of 
2 

Rp(indep) of magnitude 0,6760. The corresponding results 
2 2 

for Rp(5,indep) relative to Rp(3,indep) show that the for- 
mer averages, band-f or-band , between 6 and 8% higher than 
the latter. Comparison of FNWC and NESC independent tests, 
both at 500 mb, shows higher mean specification of the 
latter by close to 5% in favor of NESC over FNWC independent 
verification in bands 2 and 3, The NESC verification sta- 
bility seems to be dominant in bands 2 and 3 with values of 
2 ^ 

R ^ 0.85. Here again, the geographic trend previously 
noted is present. 

2 

As has been described before, R (5,indep) also has a 

r 

2 

band-2 minimum but at a level of 10% higher than R (3,indep). 

r 
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TABLE 8. Statistical results of the application of Eq . (10) to 

NESC independen t- te s t samples at 500 mb. 
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TABLE 12. Summary of NESC and FNWC independent regression test 
results. [Values are fractional explained variance (R^).] 
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Avg . 

500+300 0.7548 0.7701 0.8640 0.7962 0.7561 0.7177 0.7761 0.7499 



VI. COMPARISONS OF INDEPENDENT TO DEPENDENT RESULTS 



A summary of independent to dependent test results is 

given in Table 13. Three groups of tests are included. 

The first stratification (part a, Table 13) is averaged 

over all the bands and the data base modes. The outstanding 

aspect of part (a), column 2, is that it summarizes the 

comparative explained variances of T and T . The high 

r N 

explained variances for the dependent case at 500 mb may be 

2 

noted, as well as the relatively modest decreases in R 
(between 6 and 8%) upon testing with the independent pooled 
samples . 

Part (b) of Table 13 averages over pressure levels and 

illustrates in summary form the distinction between the NESC 

and FNWC regression samples. These samples show the geo- 

2 2 

graphic patterns of R noted earlier: e.g. R^(dep) has a 

2 

tendency to increase poleward whereas R (dep) has a band 2 

r 

minimum. The independent test column of part (b) summarizes 

the pattern of T -verification by bands (but independent of 

r 

2 

levels). Upon verification, the mean-regression R shrink- 
age of T averages between 4 and 6% when all bands are con- 
s ider ed . 

Part (c) of Table 13 summarizes some of the key ideas 

put forth in earlier sections. For example, the variable 

2 

T..(5 dep) is characterized by a high value of R (5). The 
N N 

-latter statistic decreases slightly (by about 5%) when 
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applied to independent- tes t data. The same remarks may be 

2 

applied to 1^,(5) with slightly smaller R^(5) - 0,874 noted 
in the dependent case coupled with a somewhat larger shrink- 
age (about 9%) after test upon independent data. The con- 
clusions regarding the comparative instability of the depen- 
dent test for T^(3) become clear from Table 13, part (c), 

2 

where it is seen that the values of R increase approxi- 
mately 3% from dependent to independent tests. 

The results for T_(3 dep) and T„(3 indep) are inconclu- 
sive in the form presented in Table 13(c), The whole set 
of results of part (c) indicates that little if any addi- 
tional specification accuracy is obtained from the four-day 
data base in lieu of the three-day data base. 
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TABLE 13. Summary of independent to dependent test results 



test being 
compared 


dependent ^ 
value of R 


independent 
value of 


Tp(5) 




0.8957 


0.8392 




0.8737 


0 .7798 


tJ(3) 

t;o) 




0.7298 


0 .7533 




0.7257 


0 . 7201 


T , band 
Tp, band 


1 


0.7795 


0 .7548 


1 


0. 8248 


0.7561 


T^, band 
Tp, band 


2 


0.7966 


0.7701 


2 


0.7532 


0 . 7177 


T^, band 
T^, band 


3 


0.8620 


0.8640 


3 


0 . 8212 


0 . 7761 


T^5). 4 
T^(3), 3 
T^(3), 4 
Tp(5), 3 


d ay 


0.8989 


0.8369 


day 


0.8928 


0.8369 


day 


0.7465 


0 . 7841 


day . 


0.7129 


0 .7226 


day 


0.8760 


0.7725 


Tp(5), 4 


day 


0 . 8714 


0 .7871 


Tp(3), 3 


day 


0 . 7247 


0 . 7504 


Tp(3). 4 


d ay 


0.7266 


0.6898 
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VII. CONCLUSIONS 



The results of Sections IV, V, and VI afford confirming 
evidence. For each case, high specification for occurs 
at the 500-mb level. In spite of considerable shrinkage, 

2 

there still remains a significant specification value of R 

at 300 mb. This was true also of the FNWC specifications. 

2 

The reduction in R associated with 300 mb T.. and T_ 

N F 

specifications was maximal in band 2 where more inter-di- 
urnal variability in the temperature and associated tropo- 
pause conditions occur. The NESC specifications obtained by 
Smith et al (1970) by the 16-predictor equation also fell 
off to a relative minimum at the tropopause level in their 
SIRS-A test. 

Not only is part of the observed lack of specification 
associated with the sampling of temperature extremes at or 
near the 300-mb level, but also the inter-diurnal varia- 
bility in cloud cover would affect the noise contained in 
the radiances due to the cloud-correction technique, which 
is not error-free. The errors due to cloud-correction must 
be maximal in channels 2 and 3 for which maximum dl/dln p 
weights occur in the upper troposphere. Hence, the mid- 
latitude band (band 2) which had most inter-diurnal as 
well as geographical sampling variability had the poorest 
300-rab specification. On the other hand, the 500-mb speci- 
fication and verification was adequate throughout. The 
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fact that the best 300-mb specification occurred in band 3 
(50-70N) where cloud tops are clima t ologically lower than 
in the other two bands tends to confirm this conclusion. 
Similar tests at additional levels would have provided a 
better perspective in assessing the use of the SIRS cor- 
rected radiance data in constructing vertical profiles. As 
a result, a more balanced view of the vertical resolution 
afforded by the SIRS regression techniques would have 
evolved • 

Several factors could work to improve the operational 
routine tested here: (1) to make use of a multi-level 

cloud-correction model for correction of SIRS-radiances 
(similar to that described by Smith et al 1970); (2) to 

employ a 16-predictor equation involving squares of radi- 
ances as well as linear terras; (3) to use a larger sample 
base confined to a smaller band zone, for example, 15^ 
latitude overlapping bands rather than the 20^ bands used 
here , 

It is felt that with these modifications, a SIRS 
operational unit attached to FNWC could profitably use the 
SIRS data in improving their update analysis over sparse 
data regions, and in determining the onset of significant 
circulation changes , 
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